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We study the zero temperature dynamics of Bose-Einstein 
condensates in driven high-quality optical cavities in the limit 
of large atom-field detuning. We calculate the stationary 
ground state and the spectrum of coupled atom and field 
mode excitations for standing wave cavities as well as for 
travelling wave cavities. Finite cavity response times lead 
to damping or controlled amplification of these excitations. 
Analytic solutions in the Lamb-Dicke expansion are in good 
agreement with numerical results for the full problem and 
show that oscillation frequencies and the corresponding damp- 
ing rates are qualitatively different for the two cases. 

PACS number(s): 03.75.Fi, 05.30. Jp, 32.80. Pj, 42.50.Vk 



I. INTRODUCTION 

The experimental realization of Bose-Einstein conden- 
sation (BEC) in dilute atomic gases a couple of years ago 
|J as a consequence of improved cooling and trapping 
techniques has dramatically boosted the study of ultra- 
cold atoms. Today, BEC is a widespread tool and a huge 
range of new phenomena has been investigated experi- 
mentally and theoretically, see e.g. Refs. g-f§ for recent 
overviews. In this context the interaction of BECs with 
laser light and optical lattices has been studied inten- 
sively |oj-|lO|] and effects such as the reduction of the speed 
of light by many orders of magnitude [ fTl| and the occur- 
rence of superradiance |l2| have been found. Recently the 
optical creation of vortices Jll|[l4|] has been demonstrated 
and many more intriguing effects have been theoretically 
predicted. 

In parallel, the field of cavity quantum electrodynam- 
ics, which studies the interaction of matter with one or 
a few single light modes, has reached such a level of so- 
phistication that the interaction of light with the internal 
and external degrees of freedom of a single neutral parti- 
cle can be observed and controlled in a very precise way 
iPI - For optical fields trapping and cooling of a single 
atom in a cavity mode has been demonstrated Jl6|,[l7j ■ It 
is thus a logical next step to combine these two success- 
ful techniques and study the interaction of a BEC with 
a high finesse optical cavity, that is, the strong coupling 
of a far detuned optical mode to the dynamics of a con- 
densate described by its macroscopically occupied wave 
function ]18| . Several groups have already been working 
along these lines and, for example, predicted the amplifi- 
cation of matter waves Q and the occurrence of dressed 
condensates EQl. In the extreme limit one could envis- 



age a large atomic cloud trapped and manipulated with 
a single photon. 

In this work we extent our recently proposed scheme 
for cooling one or a few atoms in high-quality optical 
cavities [^l|,^2| to the case of a BEC. This requires a 
quantum mechanical treatment of the external degrees 
of freedom and the inclusion of atom- atom interactions. 
The system under investigation is a Bose-Einstein con- 
densate interacting with the mode (the two modes) of a 
driven standing wave cavity (ring cavity) . The properties 
of such a system as a measuring device for condensates 
have been discussed previously p^j2^ |. 

We will investigate the ground state and collective ex- 
citations |24|-|27|] of the coupled condensate-light field sys- 
tem in the optical potential of the cavity. Because of the 
strong coupling of the condensate wave function to the 
cavity modes, an oscillation of the condensate also leads 
to an oscillation of the intracavity light fields. Accord- 
ingly, the oscillation frequencies of the collective excita- 
tions are shifted with respect to an external optical po- 
tential of fixed depth, i.e., the optical potential formed by 
a free space standing wave. Furthermore, for appropri- 
ately chosen parameters the dissipative dynamics of the 
cavity due to cavity decay gives rise to damping of the 
condensate excitations without incoherent spontaneous 
emission. We analyze this damping mechanism and study 
its parameter dependence by numerical solutions of the 
coupled equations of motion as well as by analytic so- 
lutions of a simplified model based on the Lamb-Dicke 
expansion. 

This paper is organized as follows. In Sec. II we discuss 
the case of a condensate interacting with the single mode 
of a standing wave cavity. After presenting the set of cou- 
pled non-linear equations of motion for the condensate 
and the light field, we discuss the numerical and analyt- 
ical solutions for the ground state and the collective ex- 
citations. Section III investigates the more complicated 
situation of a condensate coupled to the two independent 
modes of an optical ring cavity. In Sec. IV we discuss the 
influence of binary collisions between the atoms within 
the condensate on the excitation frequencies and damp- 
ing rates. Finally, we summarize our results in Sec. V. 



II. BEC IN STANDING WAVE CAVITY 

Let us first consider the case of a Bose-Einstein con- 
densate interacting with a single standing wave mode. 
The cavity mode is assumed for all times to be in a co- 
herent state \ct(t)) and the condensate is described by a 
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single wave function \i/j(t)) for all N particles, which is 
a good approximation at zero temperature. This means 
that we factorize the quantum state of the system and 
thus neglect any entanglement between the condensate 
and the cavity field which might build up in the course 
of the time evolution. This simplification is only justi- 
fied in the limit of a large photon number \a\ 2 . To avoid 
spontaneous emission, we assume a very large detuning 
of the light field from the atomic resonance. More pre- 
cisely, we assume that the cavity decay is the dominant 
incoherent process in the system and we thus require that 
k 3> NTs, where k is the cavity decay rate, N the num- 
ber of atoms in the condensate, T the spontaneous decay 
rate of the atoms, and s the atomic saturation param- 
eter. After adiabatic elimination of the atomic excited 
states we obtain the following equations of motion, 

j t a{t) = [iAc - iN{U{x)) - k] a(t) + r h (la) 



p 
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\a(t)\ 2 U(x) + Ng coU \^(x,t)\ 2 }^(x,t). (lb) 



Here A c is the detuning of the pump field from the cav- 
ity mode, U(x) = Uocos 2 (kx) is the optical potential 
formed by a single cavity photon, and g describes the ac- 
tion of the driving laser. The expectation value (U(x)) 
has to be taken with respect to the momentary wave 
function \ip(t)). Equation (lb) is the well known Gross- 
Pitaevskii equation (GPE) for a condensate in an ex- 
ternal field which in our case depends on the momen- 
tary field intensity |a(i)| 2 . The last term in the GPE 
models the interaction of atoms within the condensate 
where g co u is related to the s-wave scattering length a by 
g C oii = 4irh a/m. 

Equations (Eh are two coupled nonlinear equations de- 
scribing the dynamics of the compound system formed 
by the condensate and the cavity field p| . The most in- 
teresting effects occur for parameters where the coupling 
between these equations significantly changes the system 
behaviour. We thus impose the condition NUq > k, 
which guarantees that the presence of the condensate 
shifts the cavity frequency efficiently into or out of reso- 
nance with the driving field. At the same time the optical 
potential depth |a| 2 /7o should be large enough to provide 
at least a few bound states for the atoms. Limitations 
and the interesting parameter regimes for this model have 
been discussed in Ref. |B3|. 



A. Ground state 

In order to obtain the ground state of the compound 
condensate-cavity system, we have to find the stationary 
solution of the system of coupled nonlinear equations (|l]). 
This can be done by elimination of |a(i)| 2 in (lb) using 
(O) , and a subsequent numerical solution of the resulting 



non-linear equation for the ground state wave function 
with the method of steepest descent. This consists of 
a numerical propagation of the GPE in imaginary time 
t = it until the wave function converges to a stationary 
state. 

In this work we will concentrate on the case of Uq > 
where the potential minima coincide with the antinodes 
of the field (low-field seeking atoms). The ground state 
wave function will thus be localized at the field antin- 
odes, thereby minimizing the coupling of the condensate 
to the light field. For a cavity resonant with the driving 
field, this means that the photon number is maximum 
for the stationary ground state. Any excitation of the 
condensate will then lead to a smaller cavity field. 

As expected we find that the ground state wave func- 
tion becomes better localized for stronger driving fields 
r\ and larger optical potentials U$. On the other hand, a 
strong atom-atom repulsion (large positive values of g C oii) 
increases the width of the BEC wave function and thus 
counteracts the confining effect of the potential. This in 
due course leads to an increased coupling of the BEC to 
the cavity field and hence a smaller field intensity. A 
more detailed analysis of the ground state wave function 
ipo(x), its energy \i, and the stationary field intensity 
|ao| 2 has been given in Ref. |E3[. 



B. Collective excitations 

Let us now turn to weak excitations of the condensate 
from the ground state. First, we will calculate the spec- 
trum of collective excitations of the condensate. In con- 
trast to fixed external fields, the trapping potential in the 
cavity depends on the BEC wave function. Hence, exci- 
tations include small deviations of the wave function and 
the cavity field a from their respective stationary state. 
We may thus write ip(x, t) = exp(— ifd)[ipo(x) + Sip(x, t)} 
and a(t) = ao + Sa(t). For convenience we have already 
included the ground state time evolution into the ansatz 
for the wave function here. Inserting this into Eqs. (|l|) 
and linearizing in dip and 8a we obtain 

i^-Ja = [-A c + N{i>o\U(x)\il) ) -m]Sa 
at 

+Na {8ip\U{x)\^ ) +Na {iPo\U(x)\6i>), (2a) 

= 1 1" + l^ol 2 ^ + 2Ng coU \iP \ 2 - A 5*P 
at { 2m J 

+Ng coU 7p 2 S^* + a Uip 8a* + a* Uip 8a. (2b) 

For large k (more precisely, for 1/k much smaller than 
the time scale of the condensate motion) , the cavity field 
follows adiabatically the changes of the wave function 
and thus 8a can be adiabatically eliminated. In this case 
one recovers the limit of Ref. . 

In general the linearized time evolution couples the de- 
viations Sip and 8a also to their complex conjugates. In 
order to obtain excitation eigenstates, i.e., periodic solu- 
tions, we thus have to use the simultaneous ansatz 
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. " [e- lut 6il*+(x) +e lut 5iP-(x) 



-iut t 

- wt 5a A 



8t]) = e" 7 * [e- lut 8i]j + {x) + e" A > 
5a = e" 7 * [e~ lvt 5a + + e iut Sa*_\ . (3) 

The collective excitations are thus defined as the solu- 
tions of the eigenvalue problem 

/ 5a+ \ / 5a+ 

5a- J _ _ . [ 5a- 

w 5ip+{x) iV1 

\ 5ip_(x) J \ 5ip-(x) 



(4) 



where M is easily obtained from Eqs. (g) as a non- 
Hermitian matrix. The complex eigenvalues have the 
form u> n = v n —ijm where v n is the oscillation frequency 
of the nth collective excitation and 7„ the corresponding 
damping rate. Note that, depending on the parameters, 
negative damping rates are possible, leading to an expo- 
nential growth of the collective excitations. In this case 
the assumption of small deviations from the ground state 
imposed above only holds for very short times. Hence by 
changing some cavity parameters we can switch between 
stable and unstable cases and generate controlled exci- 
tations of the condensate and study their decay. In the 
following we will, however, concentrate on the case of 
positive 7„ and therefore damped excitations. 

Physically this damping arises from a kind of Sisyphus 
mechanism. For cavity damping rates k of the order of 
the oscillation frequencies u n , the cavity field follows with 
a certain delay the changes of the condensate wave func- 
tion. By properly choosing the system parameters, it 
can be achieved that on average the wave function has to 
climb up the potential hills at higher cavity field intensi- 
ties and runs down at lower intensities. The condensate 
thus loses potential energy which is carried away by the 
cavity output field without an intrinsic decoherence of 
the condensate. 

Furthermore it should be emphasized that the appear- 
ance of a damping rate in the linearized equations (|2|) is a 
purely quantum feature related to the width of the atomic 
wave function. In the semiclassical limit of a point-like 
particle, the self-consistent ground state yields a particle 
exactly located at the antinodes of the cavity and hence 
all expectation values in Eq. (|2a|) vanish. Thus the cavity 
field decouples from the atomic degrees of freedom and 
no damping of the atomic motion occurs to lowest order 
in the elongation x. This is in contrast to the case of a 
ring cavity as will be shown in the following section. 

In Fig. m we show the oscillation frequencies and damp- 
ing rates of the lowest collective excitations obtained nu- 
merically by calculating the eigenvalues of Eq. (Ej) on a 
spatial grid. The eigenvalues are plotted as a function 
of the cavity decay rate k. Note that in order to keep 
the optical potential constant, we also have to scale the 
driving field rj 2 and the optical potential per photon Uq 
proportional to k. 

We see that there exists one single eigenvalue uit = 
Vf —i'jf which scales approximately proportional to n in 
contrast to all of the other eigenvalues. This specific exci- 
tation of the system corresponds to an eigenmode where 



mainly the cavity field oscillates and the condensate wave 
function is only weakly perturbed. In fact, equation ( |2a| ) 
shows that in the case where the atoms are well localized 
at the antinodes of the field (semiclassical limit), the cav- 
ity mode decouples from the matter wave function and 
the eigenvalue is given by u) t — — A c — in. 
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FIG. 1. Collective condensate excitations: (a) frequencies 
v n and (b) decay rates 7„ vs cavity decay rate k. The param- 



eters are g c 



A c = 0, NU = 10k, ?f = 20kNujr. 



Second, we notice that out of the other modes v n , 
n = 1, 2, . . ., the ones with odd indices are independent 
of k and their damping rates vanish. This effect is due 
to the spatial symmetry of the problem considered here. 
For all parameters we find that the ground state wave 
function ipQ is symmetric in the position x. Thus for all 
antisymmetric excitations the expectation values in Eq. 
( pa] ) vanish and the light field decouples. Therefore these 
odd (n = 1, 3, . . .) excitations are the same as for a trap 
of constant light intensity and hence there is no Sisyphus 
damping mechanism at work. Consequently only the low- 
est symmetric collective excitations are significantly al- 
tered by the interaction with the damped cavity mode. 
We will discuss the parameter dependence of the excita- 
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tion n — 2 by using an approximate analytic solution in 
the next section. 

Let us finally emphasize that the oscillation frequency 
and the damping rate of the symmetric collective exci- 
tations can be monitored nondestructively via the cavity 
output intensity. 



C. Harmonic oscillator approximation 

In order to gain more insight into the parameter depen- 
dence of this damping mechanism, we will now analyti- 
cally solve an approximate model of our system. To this 
end we expand the optical potential U(x) = Uq cos 2 (kx) 
with Uq > up to second order around the antinodes 
of the field, i.e., we set U(x) = Uo(kx) 2 and assume 
9coii = 0, i.e., no atom-atom interaction. For simplic- 
ity we will also assume A c = 0. 



The ground state of the Schrodinger equation (lb) is 
thus the well known harmonic oscillator ground state 
which depends on the cavity field \a\ 2 in a parametric 
way. After inserting this wave function in the expec- 
tation value in Eq. (|la|) we obtain an equation for the 
self-consistent cavity field with the solution 



"0 = 9 



ir_ _ n z u lor 

2 4^.2 



(5) 



The corresponding harmonic oscillator frequency is then 



lo q = 2us R yJ\a \ 2 U /<jj R 



(6) 



and the ground state energy is /i = Wq/2. 

For the collective excitations we now have to solve Eqs. 
(||) with the harmonic potential. The last expectation 
value of Eq. (Eaf) thus reads 



{i>o\U(x)\5fl>) =U (M(kx) 2 m 
= -U —{ip \(a 

UJQ 



a^) 2 m, (7) 



where we have used the standard relation between the 
position operator x and the ladder operators a and a 1 
of the harmonic oscillator. From this we see that the 
cavity field only couples to wave function deviations Sip 
containing the ground state tpo and/or the second excited 
state ip2 of the harmonic oscillator. Most of the harmonic 
oscillator excited states are thus unperturbed and we 
find the collective excitations of the form 8ip+(x) = ijjn, 
Sa+ = <5a_ = 6ip-(x) = with the positive eigenvalues 
lu = nujQ for n = 1 and n > 3. Analogously there exist 
excitations with negative eigenvalues u — -nwo of the 
form 5ip-(x) = ip n , 5a + = 5a- = SiJj + (x) = 0. Hence, 
in addition to the antisymmetric states already found 
to decouple previously, also the higher lying symmetric 
states decouple in the harmonic approximation. There- 
fore these symmetric excitations are only damped due 
to the anharmonicity of the potential and due to atomic 
collisions in the full model. 



The remaining (and most interesting) collective exci- 
tations are finally found by restricting the wave func- 
tions Si(j± in Eq. (||) to the 2-dimensional Hilbert space 
spanned by ipo and ip2- The resulting 6x6 matrix has two 
zero eigenvalues and the other four eigenvalues have to 
be found by solving the 4th order polynomial equation 



(-iK + JVl/o^-w)(-»« 
-A{NU u R f = 0. 



u u 



lo)(u 2 - 4Wq) 



(8) 



This gives us the (complex) eigenvalues ujf and u>2 and 
their counterparts of negative frequency. Although an 
analytic solution of (||) is possible in principle, the re- 
sulting expressions are rather long and do not provide 
much insight. Instead, we calculate the eigenvalue L02 in 
the limit of large k as in Fig.|l|, i.e., by keeping Uq/k and 
rj 2 /k constant. The zeroth order in this expansion in ui/k 
yields the leading order of the frequency 



v 2 = 2u> \ 1 



N 2 U uj r 



4rf 



(9) 



and the first order gives the leading order of the decay 
rate 
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4N 2 U 2 uj 2 



1 - 



N 2 U m R 
i-q 2 



(10) 



Equation (JsJ) gives a quantitative explanation for the fre- 
quency shift of vi according to the coupling of the BEC 
and the cavity mode as compared to the value 2wo for 
the case of a harmonic oscillator potential of fixed pho- 
ton number. We also see that the small variation of z/ 2 
in Fig.|l|(a) for small values of k are in fact of the order 
1/k 2 . Equation ([To]) leads to the asymptotic behaviour 
like 1/k for the decay rate 72 in Fig.^J(b). In the limit 
of a strong driving field the frequency v 2 of the second 
collective excitation approaches the harmonic oscillator 
value. Simultaneously the damping rate 72 tends towards 
a constant non-vanishing value which is proportional to 
the square of the atom number TV. A higher condensate 
density thus significantly increases the damping of the 
collective excitation. 



III. BEC IN A RING CAVITY 

In this section we will now discuss the case of a BEC 
in a ring cavity. In this case the condensate is coupled to 
the two independent travelling wave modes a± . For sim- 
plicity we will assume in the following that both modes 
of the cavity are driven with the same pumping rate rj 
p8| . Therefore the equations of motion read 



— a±(t) = [iA c - iNUo - k] a±(t) 



(11a) 



4 



i -1p(x,t) 



A 2 

I2r7 



5^ 
2m 
+[/ |a+(t)e 



Ng coU \^{x,t)Y 

ikx i _ „—ikz 



a- (t)e 



(lib) 



Hence in general the condensate will scatter light between 
the left and right running waves and induce a strong 
coupling. This gives additional degrees of freedom to 
the system compared to the standing wave case. For ex- 
ample, in addition to intensity shifts the condensate can 
also induce a relative phase shift between the two modes, 
which changes the position of the potential wells. Anal- 
ogously, the minima can also be controlled externally by 
the relative phase of the two driving fields which allows 
to selectively excite antisymmetric excitations. 

Considering the important role which the spatial sym- 
metry plays for the standing wave cavity, we will now 
change to the description of the cavity modes by a s = 
a+ + a_ and a a — ct+ — a_ which have symmetric and 
antisymmetric mode functions, respectively. In this new 
basis Eqs. ([ll]) read 

^jOt a {t) = [iA c - iNU - iNU (cos(2kx)) - n] a s (t) 



+NU (sm{2kx))a a (t) + 2m 



—aJt) = [iA c - iNU + iNU (cos(2kx)) - k] a a (t) 
dt 



-NU (sm(2kx)}a s (t), 



(12a) 
(12b) 



i^-ip{x,t) = \^-+Ng coll \^{x,t)\ 2 
at 1 2m 

+Uq \a s (t) cos(kx) + ia a (t) sin(fca;)| 2 | 

Xip(x,t). (12c) 

Note that because of the assumption of a single pumping 
rate rj for a + and a_, in the new basis only the sym- 
metric mode a s is pumped. The antisymmetric mode a a 
only contains photons which have been scattered by the 
condensate out of a s . 

A. Ground state 



For the calculation of the ground state of the com- 
pound system formed by the BEC and the cavity modes 
we will again assume the case U > 0. We then find that 
the ground state wave function is localized at the antin- 
odes of the driven mode a s and is symmetric in x. Thus 
the expectation values of sin(2fcx) in Eqs. ( [l2"| ) vanish and 
Eq. (12b) decouples. The stationary state of the antisym- 
metric mode is therefore given by a a .o = 0. Equations 
(12a) and (12c) then reduce to the equations ([!]) for the 
standing wave cavity if one identifies the parameters 



a 



U,s|2 



2UI 
| 2 /2 
v/2V = Tf , 



for the ring cavity and the standing wave cavity, respec- 
tively. The ground state of the system can thus be ob- 
tained by using our previous results for the standing wave 
cavity and all of the discussions there equally apply to 
the ground state in the ring cavity. 



B. Collective excitations 

The collective excitations are calculated with the same 
method as in the preceding section by linearization of 
the equations of motion ( |12|) in small deviations of \ip), 
a s and a a from their stationary states \ipo), a s ,o and 0. 
Choosing the ground state wave function to be real and 
taking its symmetry into account we obtain 



i—oa s 
dt 



= f-A c 



2NU {il>o\ cos 2 {kx)\ip ) - in}5a. s 



1-nOOLa 

at 



i— dip 



dt 



(13) 



+2NU a Si o[(Sip\ cos 2 (fcr)|^ ) + c.c], (14a) 
[-A c + 2NU Q (ip \sm 2 (kx)\iP ) - in]5a a 
-NU a St0 [{6^\ sin(2fcx)|^ ) + c.c], (14b) 
+ Ko| 2 t/ O cos 2 (fci) + 2Ng coll \ 1 p \ 2 

I 2m 

-fi\5i() + Ng co ui)l5i>* 

+Uq cos 2 (kx)ipo( a s,o8ct* s + c.c.) 
-iU sm(2kx)Tpo(a StO 6a* a -c.c). (14c) 



From these equations we see that the behaviour of the 
excitation eigenstates strongly depends on their spatial 
symmetry. 

For sy mmet ric excitations Sip(x) the last expectation 
value in (14b) vanishes and the antisymmetric cavity 
mode decouples from the wave function. Hence in this 
case we find 5a a = 0. The equations of motion for Sip 
and 6a s then reduce to their standing wave counterpart 
discussed in the previous section if one rescales the pa- 
rameters as in (|l3|). The symmetric collective excitations 
are thus the same as those in a standing wave cavity. 

Analogously, for antisymmetric excitations Sip(x) the 
symmetric cavity mode decouples and therefore Sa s = 0. 
We then find a new set of coupled equations for Sip and 
Sa a . Thus, in contrast to the case of a standing wave 
cavity, also the antisymmetric excitations are damped 
in a ring cavity. However, the damping mechanism is 
of completely different physical origin. Instead of the 
Sisyphus mechanism discussed above, here the coherent 
scattering of photons from the a s cavity mode into the 
a a mode is responsible for the damping. This leads to 
less severe requirements for the cavity parameters as we 
will see in the following subsection. 

Figure ^ shows the spectrum of collective excitations 
of a Bose condensate in a ring cavity which is obtained 
from the numerical solution of Eqs. ([l4]). First, we note 
that in contrast to the case of the standing wave cavity 
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we now find two modes with eigenvalues which scale pro- 
portional to the cavity decay rate k. In the semiclassical 
limit (atoms well localized) , these correspond to pure os- 
cillations of the symmetric and antisymmetric field mode, 
respectively, and are thus labelled u s = v s — ij s and 
u a = v a — i"f a . The semiclassical limits of these eigcnfrc- 
quencies are obtained from Eqs. (Q) as u) s = — A c — in 
and ui a = — A c + 2iV£/o — in. Although the damping rates 
of these modes are equal, we see that the different spa- 
tially depending coupling to the atoms leads to a large 
difference in the oscillation frequencies. 
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FIG. 2. Collective condensate excitations in a ring cavity: 
(a) frequencies v n and (b) decay rates y„ vs cavity decay 
rate k. The parameters are g co ii = A c = 0, NUo = 5k, 
rf = 70kNur. 

For the parameters chosen in Fig. || the other oscilla- 
tion frequencies v n , n > 1, are mainly independent of k. 
However, whereas all frequencies with n > 2 are equally 
spaced and hence very well described by harmonic oscilla- 
tor states, the lowest frequency v\ is significantly shifted 
downwards. For the damping rates we find that only the 
two lowest excitations exhibit relevant damping. How- 
ever, the dependence of these damping rates on the sys- 
tem parameters is qualitatively different according to the 



different damping mechanisms. We will return to the dis- 
cussion of these features in the following subsection where 
we derive analytic approximations for the eigenvalues. 



C. Harmonic oscillator approximation 

Let us now calculate analytic estimates for the lowest 
oscillat ion f requencies and damping rates along the lines 
of Sec. IIC. We will thus again assume A c = g co u = 0. 

As we have already seen, the calculation of the ground 
state and of the symmetric collective excitations can be 
reduced to the problem of the standing wave if the ap- 
propriate identification of the system parameters ( |l3| ) is 
made. We can therefore use our previous results to ob- 
tain the self-consistent cavity field 



\ a s,o\ = 



Arf N 2 U u R 



and the corresponding harmonic oscillator frequency 



(15) 



uj = 2uj R ^J\a sfi \ 2 U /u) R . 



(16) 



For the lowest symmetric excitation, the expansion for 
large values of k yields 



v-2 = 2u \ 1 



N 2 U u; R 



and 
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16N 2 U$lu 



4r) 2 



1 N 2 U u R Y 
4r] 2 J 



(17) 



(18) 



Analogously, we can calculate the lowest antisymmet- 
ric excitation by expanding the expectation values in 
Eqs. (14b) and (14c) to lowest order in kx. To this or- 
der only Sip and Sip* proportional to the first harmonic 
oscillator wave function tpi couple to the cavity field Sa a 
and 5a*. and we thus have to find the eigenvalues of a 4x4 
matrix, that is, we must solve the characteristic polyno- 
mial 



-««+ 2AT/ (1 



UJO 



)_ w]Hk _2AT/ (1-^)- w ] 



x (uj 2 - u 2 ) - 4(iVcW ) 2 (l - = 0. 



(19) 



In the limit of k — * oo (with constant Uq/k and rf / ' k) 
this yields the oscillation frequency 



V\ = uo\ 1 - 



4N 2 U 2 (l-^) 

K 2 +4N 2 U 2 {1 _^ ) 2- 



(20) 



The first order correction in 1 /k gives the dominant term 
of the corresponding damping rate 
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AN 2 'V 2 (1 



Wo ' 



'[ K 2 +4 A2[/ 2 (l-^)2]2- 



(21) 



G 



We can now compare the behaviour of the two lowest 
eigenvalues as a function of the system parameters. As 
an example, let us consider the case of a relatively strong 
pump, i] 2 ^S> N 2 Uouir. In this limit, the second excita- 
tion frequency is only weakly shifted from the har- 
monic oscillator frequency 2ljq. On the other hand, the 
frequency shift of the lowest excitation v\ mainly de- 
pends on the ratio NUq/k. Since this ratio has to be 
larger than one in order to yield a significant frequency 
shift of the cavity by the atoms, equation (|2(]) implies 
that v\ is strongly shifted towards zero. Simultaneously 
we find for the damping rates that 72 becomes indepen- 
dent of r\ in this limit, in contrast to 71 which is pro- 
portional to ujq and thus proportional to i] 2 . Therefore 
the damping rate of the first antisymmetric excitation 
can be increased arbitrarily by increasing the intensity 
of the pump field. The damping rate of the first sym- 
metric excitation is much harder to manipulate because 
it is mainly governed by the optical potential per photon 
and thus by the quality of the cavity. On the other hand, 
we note that 72 scales proportional to TV 2 whereas 71 is 
inversely proportional to N 2 . The number of atoms thus 
provides another handle to change the relative size of the 
damping rates 71 and 72. 

Another point is worth a comment here. We empha- 
sized in the previous section that the damping mechanism 
for the collective excitations in a standing wave cavity is 
crucially related to the width of the matter wave function 
and vanishes in the semiclassical limit where the atoms 
are treated as point particles. In contrast to this we find 
that in the travelling wave cavity the damping mecha- 
nism still exists in the semiclassical limit. In fact, our 
results for the oscillation frequency (|2(]) and the damp- 
ing rate (21) agree with the semiclassical results 29 if 
one takes formally the limit u>r/u>q — ► 0. 

In Fig. |3| we show the excitation frequencies and 
the damping rates 71.2 as a function of the pump strength 
r} 2 for both the numerical solution and the analytic ap- 
proximations. We see that for the chosen parameters 
the approximations fit quite well apart from the values 
of 72. This comes from the fact that we obtained the 
complex eigenvalues u) n from an expansion of Eqs. (|^) 
and (|ig| ) for small values of <C 1. As we see from 

Fig. Wa) this is well fulfilled for lji for the chosen pa- 
rameters, but \u>2\/k is of the order of one. However, 
in the limit of a strong pump the lowest order term for 
the frequency already gives the correct value, namely 
twice the harmonic oscillator frequency. Thus, only the 
imaginary part (the damping rate 72) of the analytic ap- 
proximation deviates from the exact solution in Fig. |3|. 
In parameter regions where |w2|/k *C 1 we find a much 
better agreement of the two solutions. 
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FIG. 3. Collective condensate excitations in a ring cavity: 
(a) frequencies v„ and (b) decay rates 7„ vs pump strength rf . 
The parameters are g co u = A c = 0, NUo = 5k, k = 400ujr. 
The solid curves are obtained numerically, the dashed curves 
are the analytical solutions obtained in the harmonic oscilla- 
tor approximation. 



IV. INTERACTING BOSE GAS 

In the discussion so far we have omitted the effects 
of atomic interactions, as described by the collision rate 
g co u in the GPE, on the energies and damping rates of 
the collective excitations. Neglecting this has allowed us 
to obtain analytical expressions and therefore to discuss 
the parameter dependence of our results explicitly. How- 
ever, atomic collisions are known to play a crucial role in 
experimental realizations of Bose-Einstein condensates. 
We will now discuss the changes of the collective excita- 
tions according to collisions in a numerical example of a 
condensate in a ring cavity. 

We show in Fig. |j the excitation frequencies v n and the 
corresponding damping rates j n as a function of the colli- 
sion rate g co u with all other parameters fixed. The main 
effect on the stationary ground state wave function |23| ] 
of a repulsive interaction between the condensed atoms 
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is to increase the width of the wave function. Since this 
larger width also changes the coupling to the cavity held, 
we find that the steady state photon number decreases 
with increasing collision rate. Consequently, the optical 
potential becomes more shallow and the excitation fre- 
quencies decrease. However, as we can see from Fig. |] 
this argument thus not hold for the lowest (antisymmet- 
ric) excitation. Here the atomic collisions counteract the 
strong frequency shift which we found in the previous 
section and v\ slightly increases with g co u. Above a cer- 
tain threshold value for g co u the atom-atom repulsion 
gets stronger than the confining effect of the optical po- 
tential. In this case the ground state wave function is no 
longer localized and the spectrum of excitations changes 
into that of unbound particles where each excitation fre- 
quency is doubly degenerate. 



increasing g co u, 72 increases significantly over a broad 
range of values of g co u . This is related to the fact that the 
damping of the symmetric excitation 72 depends crucially 
on the width of the ground state wave function whereas 
the damping mechanism of the antisymmetric excitation 
does not, as already emphasized before. Since the major 
effect of atom collisions is to broaden the wave function, 
the resulting changes of the damping rates occur predom- 
inantly for the symmetric excitations. Note also that the 
atomic collisions and the stronger anharmonicity of the 
potential according to the lower field intensity enhance 
the damping of higher collective excitations, as can be 
seen from the damping rates 73 and 74 in Fig. ^(b) . 

V. CONCLUSIONS 
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FIG. 4. (a) Frequencies v n and (b) decay rates 7n of the 
collective excitations vs atomic collision rate g co u- The pa- 
rameters are g coU = A c = 0, NUo = 2k, k = 400lor, 
rf = 200Nklj r . 

Figure ^(b) shows that collisional effects have an even 
more important influence on the damping rates of the col- 
lective excitations. We see that the effect differs for the 
damping rates 71 and 72. While 71 weakly decreases with 



In summary we have studied in detail the interaction 
of a Bose-Einstein condensate with one or two single 
modes in a high-finesse optical cavity. We have solved 
the coupled set of non-linear equations of motion for the 
joint dynamics of the condensate and the light field nu- 
merically and compared it analytically with a simplified 
model based on the Lamb-Dicke expansion. We find that, 
even without atom-atom interaction, the oscillation fre- 
quencies are shifted with respect to their values in a fixed 
external potential. 

For a finite cavity response time the collective excita- 
tions are damped or amplified depending on the cavity 
detuning, which can be easily controlled externally. We 
identify two distinct mechanisms depending on the spa- 
tial symmetry of the excitations. The damping mech- 
anism in standing wave cavities and for spatially sym- 
metric excitations in ring cavities is due to a Sisyphus 
type effect, which leads to larger cavity fields at times 
when the condensate runs up potential hills than at times 
when the condensate runs down. On average this effect 
extracts kinetic energy from the condensate which is car- 
ried away by the cavity field. On the other hand, the 
damping mechanism for the spatially antisymmetric ex- 
citations is only present in a ring cavity due to the scat- 
tering of cavity photons between the two counterpropa- 
gating waves. This creates an intensity imbalance, which 
is counteracted by the cavity damping and hence leads 
to momentum dissipation. 

The two damping mechanisms exhibit very distinct 
parameter dependences. Our analytical approximations 
show that in the limit of strong cavity pumping the 
damping rate of the spatially symmetric excitations be- 
comes independent of the pump but scales proportional 
to the square of the atom number N 2 while the damp- 
ing rate of the antisymmetric excitations is proportional 
to the pump field intensity and inversely proportional to 
./V 2 , which implies less stringent requirements to cavity 
technology. 

The difference between the damping/amplification 
rates of excitations with different spatial symmetry could 
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be used to manipulate a Bose-Einstein condensate in a 
controlled fashion. In addition, in a ring cavity setup we 
can also excite oscillations by external phase and ampli- 
tude shifts of the pump light. All the effects could of 
course be enhanced by tailored feedback of the measured 
transmitted intensity onto the pump. This might give 
rise to useful applications of such a system in the con- 
text of quantum information and quantum computation 
in analogy to other recently proposed systems making 
use of particles in optical lattices 
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